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the  fundamental  mode.  The  Falkner-Skan-Cooke  flow,  modulated  by  the  presence  of  the 
crossflow  vortex,  is  found  to  be  subject  to  a  new  secondary  instability  with  large  growth 
rates. 


Nonlinear  Development  of  Gortler  and  Crossflow  Vortices  and 
Gortler/Tollmien-Schlichting  Wave  Interaction 

ABSTRACT 

The  problem  of  nonlinear  development  of  Gortler  vortices  on  a  curved 
wall  is  studied  within  the  framework  of  incompressible  Navier-Stokes 
equations  which  are  solved  by  a  Fourier-Chebyshev  spectral  method.  The 
results  show  that  higher  harmonics  grow  due  to  nonlinear  effects; 
however,  most  of  the  energy  remains  in  the  fundamental  mode.  The 
computed  flow  field  in  the  presence  of  a  Gortler  vortex  is  in  qualitative 
agreement  with  the  experimental  data.  The  interaction  of  the  Gortler 
vortex  with  a  two-dimensional  Tollmien-Schlichting  wave  is  also  studied 
and  it  is  shown  that  the  Tollmien-Schlichting  wave  grows  faster  than  its 
linear  theory  growth  rate  when  the  amplitude  of  the  Gortler  vortex  is  suffi¬ 
ciently  large.  Due  to  nonlinear  effects  this  interaction  further  leads  to  the 
development  of  oblique  waves  with  spanwise  wavelength  equal  to  the 
Gortler  vortex  wavelength.  The  numerical  method  is  also  applied  to  study 
the  nonlinear  development  of  a  stationary  crossflow  vortex  in  a  Falkner- 
Skan-Cooke  boundary  layer.  The  crossflow  vortex  develops  in  a  manner 
similar  to  that  found  earlier  for  rotating  disk  flow.  The  fundamental  and 
the  higher  harmonics  all  tend  to  saturate  when  the  integration  is  carried  to 
large  amplitudes.  The  computed  velocity  distribution  clearly  shows  the 
emergence  of  the  superharmonic  which,  however,  does  not  dominate  the 
fundamental  mode.  The  Falkner-Skan-Cooke  flow,  modulated  by  the 
presence  of  the  crossflow  vortex,  is  found  to  be  subject  to  a  new  secondary 
instability  with  large  growth  rates. 
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1. 


The  subject  of  laminar/turbulent  boundary  layer  transition  is  of 
fundamental  and  practical  importance  in  fluid  mechanics.  Depending 
upon  the  flow  conditions  and  body  geometry,  various  instability  mecha¬ 
nisms  such  as  Tollmien-Schlichting  (TS),  Gortler  and  crossflow  may  be 
operative  and  could  also  interact  if  multiple  instability  modes  are  present. 
A  number  of  studies  have  been  carried  out  to  understand  the  breakdown  of 
TS  waves  to  turbulence  in  a  flat-plate  boundary  layer.  Notable  among  these 
are  the  experiments  of  Klebanoff  et  al.  [1],  Kachanov  et  al.  [2],  Saric  & 
Thomas  [3],  the  theoretical  studies  of  Craik  [4]  and  Herbert  [5],  and  numer¬ 
ical  simulations  of  Wray  &  Hussaini  [6],  Orszag  &  Patera  [7],  Spalart  & 
Yang  [8],  Laurien  &  Kleiser  [9],  Zang  &  Hussaini  [10]  and  Fasel  et  al.  [11]. 
Because  of  these  and  other  related  studies,  at  least  the  early  stages  of  the 
breakdown  process  are  now  relatively  well-understood. 

The  boundary  layer  that  develops  on  a  concavely  curved  surface  is 
subject  to  Gortler  instability  which  manifests  itself  in  the  form  of  steady 
streamwise  counter-rotating  vortices.  Similarly,  near  the  leading  edge  of  a 
swept  wing,  crossflow  instability  sets  in  due  to  inflectional  crossflow 
velocity  profiles.  This  instability  is  known  to  result  in  corotating  steady 
vortices  known  as  crossflow  vortices  which  are  oriented  at  small  angles 
with  the  inviscid  streamlines.  Sometimes,  this  instability  may  also  result 
in  unsteady  crossflow  vortices.  Depending  upon  flow  parameters, 
Gortler/TS  and  crossflow/TS  modes  may  coexist  and  mutually  interact. 
This  problem  has  received  very  little  attention,  even  though  it  is  of  direct 
application  in  the  design  of  modem  aircrafts. 

Experimental  work  on  Gcrt’cr  vortices  has  been  performed  for 
example,  by  Bippes  &  Gortler  [12],  Winoto  &  Crane  [13]  and  Swearingen  & 
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Blackwelder  [14].  The  work  of  the  latter  authors  was  particularly  carried 
out  to  understand  the  details  of  the  growth  and  breakdown  of  Gortler 
vortices  and  eventual  transition  to  turbulence.  The  results  showed  the 
development  of  low  speed  regions,  between  the  counter-rotating  vortices, 
which  project  outside  the  nominal  Blasius  boundary  layer  and  create 
strong  inflectional  normal  and  spanwise  velocity  profiles.  These  profiles 
develop  secondary  instabilities  resulting  in  unsteady  oscillations  prior  to 
transition.  Two  types  of  secondary  instability  were  identified.  The  first 
results  in  a  horse-shoe  vortex  structure  also  found  by  Aihara  &  Koyama 
[15].  The  second  type  of  instability  results  in  sinuous  transverse  oscillations 
which  were  earlier  reported  by  Bippes  [16].  Both  modes  were  also  observed 
in  the  experiment  of  Jeans  &  Johnston  [17] 

Grossflow  instability  in  swept-wing  flows  was  first  discovered  by 
Gray  [18]  and  later  the  work  of  Owen  &  Randall  [19]  and  Gregory  et  al.  [20] 
made  the  phenomenon  more  clear  and  also  provided  a  theoretical  basis  for 
appearance  of  this  instability  in  general  three-dimensional  boundary 
layers.  Experimental  investigations  into  this  instability  in  rotating  disk 
flow  were  performed  by  Smith  [21],  Gregory  et  al.  [20]  and  Wilkinson  & 
Malik  [22],  among  others.  Due  to  three-dimensional  basic  flow,  questions 
still  remain  about  the  nature  of  the  instability  even  in  the  linear  regime. 
Historically,  it  had  been  believed  that  this  instability  would  only  result  in 
steady  crossflow  vortices.  However,  the  work  of  Malik  &  Poll  [23]  pointed 
out  that  small  unsteady  disturbances  may  also  play  an  important  role  in 
the  transition  process  in  these  flows.  Since  then,  a  number  of  studies  have 
been  initiated  to  understand  the  nature  of  the  swept-wing  boundary  layer 
instability  at  its  linear  and  nonlinear  stage.  Notable  among  these  are  the 
experiments  of  Muller  &  Bippes  [24]  and  Dagenhart  et  al.  [25].  Both  the 
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stationary  and  unsteady  disturbances  were  observed  in  these  experiments 
which  concern  low-speed  flows.  Earlier,  Saric  &  Yeates  [26]  only  found 
stationary  crossflow  vortices.  However,  they  observed  that  at  some  distance 
downstream  from  the  leading  edge  a  superharmonic  appears  in  the  hot¬ 
wire  signal.  The  amplitude  of  the  superharmonic  was  found  to  be  three 
times  the  amplitude  of  the  fundamental.  Reed  [27]  attempted  to  explain  this 
phenomenon  by  using  her  wave-interaction  theory.  Numerical  simulations 
of  the  early  nonlinear  stages  of  the  crossflow  vortex  have  been  performed  by 
Malik  [28]  and  recently  by  Meyer  &  Kleiser  [29]. 

The  Gortler/TS  interaction  problem  has  been  studied  by  Malik  [30] 
using  numerical  methods.  Hall  &  Smith  [31]  using  asymptotic  methods  and 
Nayfeh  &  Al-Maaitah  [32]  using  the  method  of  multiple  scales.  Questions 
still  remain  that  at  what  amplitudes  such  interactions  take  place.  An 
understanding  of  these  interactions  will  also  shed  some  light  on  the 
crossflow/TS  interaction  problem.  The  practical  interest  in  these  interac¬ 
tion  studies  arises  from  their  relevance  in  boundary  layer  transition  predic¬ 
tion  which  is  needed  for  the  design  of  modem  aircrafts.  The  state  of  the  art 
technique  used  for  transition  prediction  is  based  upon  linear  theory  in  the 
form  of  the  eN  method.  If  wave-interactions  of  the  type  discussed  above  were 
to  take  place,  they  will  render  the  prediction  of  the  eN  method  useless  and 
will  have  disastrous  consequences  for  the  design  of  laminar  flow  control 
wings.  It  is,  therefore,  important  to  understand  the  nature  of  these  wave 
interactions  and  determine  the  parameter  space  within  which  they  might 
take  place. 

In  the  present  report,  we  use  a  Fourier/Chebyshev  spectral  method  to 
investigate  the  nonlinear  development  of  a  Gortler  vortex.  The  computed 
flow  field  is  compared  with  the  experiments  of  Swearingen  &  Blackwelder 
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[14].  The  numerical  method  is  then  used  to  study  the  interaction  of  a  small 
amplitude  TS  wave  with  a  finite-amplitude  Gortler  vortex.  Finally, 
nonlinear  development  of  a  steady  crossflow  vortex  is  considered  and  the 
emergence  of  a  secondary  instability  is  discussed. 


2.  GOVERNING  EQUATIONS  AND  NUMERICAL  APPROACH 

We  consider  the  problem  of  laminar/turbulent  transition  in  two-  and 

three-dimensional  boundary  layers.  The  governing  Navier-Stokes 

equations  for  an  incompressible  fluid  are: 

Ut+U-VU  =  -Vp+  vV2U  +  F,  (1) 

V-t/  =  0,  (2) 

where  U  =  (u ,  v ,  w)  is  the  velocity  vector,  p  the  pressure,  and  V  the 
kinematic  viscosity.  The  term  F  is  a  forcing  term  to  be  specified  later.  Slip- 

free  boundary  conditions  at  a  solid  wall  (y  =  0)  are  imposed,  i.e., 

U  ={x,0,z,t)  =  0. 

In  the  free  stream,  it  is  reasonable  to  require  that 

U-*U„  as  y  -*  °o.  (4) 

Equations  (1-4)  are  solved  by  a  Fourier-Chebvshev  spectral  method. 
Periodicity  is  assumed  in  the  x  and  z  directions.  The  dependent  variables 
have  Fourier-Chebyshev  series  of  the  form 


u(x,y,z,t)=  ^  £  £u(Oe2,o(M//')e2”(v/4)rft 


^  2  ^  2 


(n) 


(5) 


where  Lx  and  Lz  are  the  periodicity  lengths  in  x  and  z  directions, 
respectively,  and  Tn  are  the  Chebyshev  polynomials  of  degree  n.  The 
spatial  discretization  employs  spectral  collocation.  The  collocation  points 
for  the  periodic  directions  are 
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Xj  =  jLx  IKX,  j  =  o,  1 . Kx-  1  (6) 

zm=mLz/Kz,  m  =  Q,l,...,Kz-l.  (7) 


A  staggered  grid  is  employed  in  the  normal  direction.  Velocities  are  defined 
at  the  points 


Vn  =  cos 


f  Ttq  ^ 


q  =  0,X...,N 


(8) 


\Nj 

and  the  pressures  at 

Vq+i/2  =  cos(/r(q+l/2)/lV),  q  =  0,1,..., IV  -  1  (9) 

No  artificial  pressure  boundary  conditions  are  therefore  needed.  The 
momentum  equations  are  imposed  at  the  points  given  by  Eq.  (8)  and  the 
continuity  at  those  given  by  Eq.  (9). 

The  normal  computational  coordinate  T)  is  related  to  y  through  the 

algebraic  transformation 
1+  T] 


y  =  a 


,  2a 
1  + - rj 

3L 


(10) 


where  a  is  a  scaling  constant  used  for  proper  distribution  of  points  and  y„  is 
the  location  where  free  stream  asymptotic  boundary  conditions  [33]  are 
imposed. 

In  the  spectral  collocation  method,  spatial  derivatives  of  u  are 
obtained  by  differentiating  the  series  expansion  coefficients  u(t)  determined 
by  discrete  Fourier  and  Chebyshev  transform  of  the  grid-point  values  of  u. 
The  temporal  discretization  involves  Crank-Nicolson  on  the  pressure 
gradient  and  vertical  diffusion  terms.  The  remaining  terms  in  the  momen¬ 
tum  equations  are  handled  explicitly  by  using  second-order  Adams- 
Bashforth  method.  The  incompressibility  constraint  is  imposed  implicitly. 
The  resulting  implicit  equations  are  solved  iteratively  by  a  minimum 
residual  method.  The  details  of  the  procedure  are  given  by  Malik  et  al.  [33]  . 
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2.1  Flow  over  a  curved  wall 

The  numerical  procedure  outlined  above  is  applicable  to  the  study  of 
transition  problem  in  two-  and  three-dimensional  boundary  layers  where 
TS,  Gortler  or  crossflow  instability  mechanisms  may  be  operative.  In  order 
to  explain  the  numerical  approach,  we  describe  the  application  of  the 
method  to  the  flow  over  a  concavely  curved  plate  where  both  Gortler  and  TS 
waves  might  exist.  Extension  to  three-dimensional  boundary  layers  will  be 
considered  in  a  later  section. 

Consider  incompressible  flow  with  free  stream  velocity  U„  along  a 
mildly  curved  wall  with  constant  curvature  K=  l  /  r  where  l  is  a  character¬ 
istic  length  scale  and  r  is  the  radius  of  curvature  of  the  wall.  If  x  is  the 
distance  along  the  curved  wall,  z  along  the  span,  and  y  normal  to  the  wall, 
then  the  governing  equations  represented  in  body-oriented  coordinates  are 


du  .  du  du  du  .dp  i 

+  Au~ -  +  v—  +  w—  +  yuu  =  -A-^-  +  v 
dt  dx  dv  dz  dx 


V2u-y2u  +  2yA  — 
dx 


(11) 


du 

dt 


du  du 


du 


dp 


+  Au~  +  u-^-  +  w-^--  yu2  -  ~-^-+  vi  V2u-  y2u  -2yX~ 


dx  dy 


dz 


dy 


du 


dx 


dw  .  dw  dw  dw  dp  „•> 

—  +  Au  —  +  u  —  +  w—  =  -—- f  vV2<y, 
dt  dx  dy  dz  dz 


.  du  du  dw 
dx  dy  dz 


where 


2  ■  to  +  r~ 


dy 


A  =  ■ 


1 


1+  Ky 


Y  = 


1+  Ky 


(12) 


(13) 


(14) 
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In  Eqs.  (11-14),  if  all  velocities  are  assumed  to  be  scaled  by  t/_,  lengths 
by  i,  time  by  l/U„,  pressure  by  pUt ,  then  V  =  1/R,  where  R  is  the  Reynolds 
number.  A  Gortler  number  can  be  defined  as  G  = 

According  to  first-order  boundary  layer  theory,  the  above  equations 
reduce  to  the  Blasius  flow  which  provides  the  basic  flow  for  our  study  of 
nonlinear  development  of  Gortler  vortices  and  Gortler/TS  interaction  in  a 
two-dimensional  boundary  layer. 

2.2  Linear  Stability  Theory 

We  consider  an  infinitesimally  small  disturbance  superposed  on  the 
Blasius  flow.  For  instance,  the  streamwise  velocity  u(x,y,z,t)  may  be 
written  as 

u(x,y,z,t)  =  UB(y)  +  e^i(y)e‘(“r+^"“<,  (15) 

under  the  quasi-parallel  flow  assumption.  Here  e  is  a  small  parameter,  a 
and  p  are  disturbance  wave  numbers  in  x,z  directions  respectively  and 
(0  =  wr  +  io)l  is  the  complex  frequency.  The  real  part  o)r  is  the  actual  distur¬ 
bance  frequency  and  the  imaginary  part  <y,  is  the  temporal  growth  rate. 
The  disturbance  amplitude  increases  if  w,  >  0,  otherwise  the  disturbance 
dies  out. 

Linear  stability  equations  may  be  derived  by  substituting  u(x,y^z,t) 
from  Equation  (15)  and  similar  expressions  for  other  dependent  variables 
into  the  Navier-Stokes  equations  (11-14)  and  then  retaining  only  terms  of 
0(e).  The  resulting  sixth-order  system  of  ordinary  differential  equations 
describes  an  eigenvalue  problem  for  parameters  a,  P  and  co  which  is  solved 
by  the  fourth-order  accurate  compact  difference  scheme  of  Malik  et  al.  [34], 
Depending  upon  the  values  of  Gortler  number  G  and  wave  numbers  (a, p), 
the  governing  equations  yield  solution  for  both  the  Gortler  vortices  and  TS 
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waves.  These  solutions  are  needed  to  provide  initial  conditions  for  the  non-  j 

linear  calculations.  In  figures  (1-2),  we  provide  linear  results  for  Gortler 
and  TS  waves  to  emphasize  differences  between  the  two  instabilities. 

Description  for  figures  (1-2)  is  provided  below. 


2.3  Nonlinear  Simulations 

Various  nonlinear  problems  may  now  be  formulated.  Suppose,  we 
are  interested  to  studv  the  nonlinear  development  of  Gortler  vortices.  A 
steady  Gortler  vortex  is  characterized  by  cor  -  ar  =  0,  G  >  0,  ft  >  0.  We  solve 
linear  stability  problem  for  Blasius  boundary  layer  and  find  f3  =  (30,  the 
growth  rate  0)t  and  the  disturbance  eigenfunctions  uG,vG,wG.  These  are 
presented  in  Figure  1.  The  plots  show  that  the  streamwise  component  uG  is 
an  order  of  magnitude  larger  than  the  cross-stream  components  uG ,  and 
ii'a.  The  spanwise  component  is  90°  out  of  phase  with  uG  and  vG .  This 
allows  steady  Gortler  vortex  problem  to  be  represented  as  a  system  of  real 
equations.  In  the  present  case,  our  complex  system  of  equations  automati¬ 
cally  yields  the  real  solution  for  the  steady  Gortler  vortex.  The  structure  ol 
uG  and  wG  is  such  that  it  results  in  a  pair  of  counter-rotating  vortices.  Here 
we  only  consider  steady  Gortler  vortices  since  in  the  experiments  these 
vortices  are  commonly  found  to  be  steady.  The  Gortler  vortex  is  assigned  an 
initial  finite-amplitude  ec.  Then  an  initial  velocity,  say  for  .r  component  of 
velocity,  is  of  the  form 

u(x,y,z,0)  =  U„(y)-f  ^Real^ty)?1^] 

Equations  (11-14)  are  solved  with  these  initial  conditions  and  solution 
is  marched  in  time.  The  solution  yields  the  development  of  the 
fundamental  and  various  harmonics  depending  upon  the  number  of  modes 
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allowed  in  the  simulation  according  to  Eq.  (5).  The  temporal  development  is 
related  to  the  spatial  development  through  group  velocity  transformation 
(Cg  =  da> /  da)  as  for  TS  waves. 

The  eigenfunctions  for  a  TS  wave  are  presented  in  Figure  2.  Th^ 
eigenfunctions  are  now  complex  with  disturbance  peak  lying  closer  to  the 
solid  boundary.  The  interaction  of  a  Gortler  vortex  with  two  oblique  TS 
waves  may  be  sludied  by  choosing  the  following  initial  conditions: 

u(x,y,z,0 1  =  Ur(y)  +  eGReal[uG(y)e!,v]  +  ers  jRealJu(y)e‘(“reX+/w)J 

+Real[u(y)e‘(arar'^)]}  ^ 

Here  various  values  of  fiTS  in  relation  to  may  be  assigned.  For 
example,  prs  =  (jG,2  describes  a  subharmonic  resonance  and  pTS  =  fiGS 
describe:  a  fundamental  resonance.  Relative  magnitudes  of  e0  and  eTS 
describes  various  regimes  of  interaction  where  either  Gortler  vortex  or  TS 
wave  is  small  or  when  both  have  large  amplitude  and  interact  nonlir, early. 

How  the  numerical  simulation  may  be  used  to  study  the  intermodal 
interaction  may  be  best  described  by  providing  the  example  of  snbharmonic 
instability  in  a  flat  plate  boundary  layer.  This  instability  which  is  the  result 
of  parametric  resonance  between  a  two-dimensional  TS  wave  and  oblique 
TS  waves  with  streamwise  wavelength  twice  that  of  the  fundamental  was 
first  theoretically  studied  by  Herbert  [5]  using  Floquet  analysis.  Here  we 
consider  fhe  fundamental  wave  to  have  a  =  .2033  and  prescribe  R  =  606. 
amplitude  of  .01  (i.e.,  1%)  is  assigned  to  the  two-dimensional  TS  wave. 

According  to  Herbert,  a  strong  subharmonic  instability  of  wave 
numbers  ((a/2),  f5)  should  develop  for  a  wide  range  of  /?,  (see  F'^.  9  in  [5]). 
The  present  calculation  is  performed  by  including  in  the  initial  conditions 
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two  oblique  modes  (a/2,  ft  =  ±  .16968).  Results  of  the  calculation  are 
presented  in  Fig.  3.  In  this  figure,  the  natural  logarithm  of  the  energy 
history  of  the  modes  (a,0)  and  (a/2,  (5)  are  plotted  as  a  function  of  time. 
Calculations  are  carried  out  to  about  3.5  linear  time  periods  of  the  finite- 
amplitude  two-dimensional  wave.  The  solid  curves  are  the  linear  theory 
results  for  the  corresponding  modes.  The  finite-amplitude  two-dimen¬ 
sional  wave  follows  its  linear  growth  (coi  =  .7E-5)  curve  within  the  computa¬ 
tional  domain.  The  subharmonic  (a/2,  p)  follows  the  linear  theory  result  for 
primary  subharmonic  for  some  time  and  then  diverges  from  it  indicating 
strong  secondary  instability.  The  growth  rate  at  the  onset  of  this  instability 
is  (o,  ~  .0098  which  drops  to  CO,  ~  .0071  towards  the  end  of  the  computation 
due  to  nonlinear  effects.  This  is  in  good  agreement  with  the  prediction  of 
Herbert  for  co,  =  .00824.  The  subharmonic  secondary  instability  has  been 
linked  to  the  squire  mode  by  Herbert.  However,  in  order  to  capture  the 
instability  using  numerical  simulations,  we  need  not  provide  the  actual 
eigenfunctions  for  the  squire  mode  in  the  initial  conditions.  In  our  calcula¬ 
tions  we  have  used  the  eigenfunctions  for  the  primary  subharmonic.  In 
another  calculation  the  initial  conditions  for  the  (a/2 ,±p)  modes  were 
arbitrarily  taken  to  be 

u  =  0,  i>  =  y2e"Avl)  ,  and  w  =  ±  —  — 

ns) 

The  results  are  presented  in  Figure  4.  Again  the  subharmonic  secondary 
instability  emerges.  Most  interestingly,  both  calculations  converge  towards 
the  same  growth  rate  for  the  subharmonic.  Higher  harmonics  also  grow  in 
these  nonlinear  calculations.  In  a  transition  simulation,  one  would  like  to 
assign  a  small  random  distribution  of  energy  to  all  the  modes,  except,  of 
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course,  the  finite-amplitude  disturbance  and  let  the  dominant  pattern 
evolve  by  itself. 

In  our  calculations  we  allow  the  mean  flow  to  develop  in  time  in 
order  to  simulate  the  development  of  Blasius  profile  along  the  curved  plate. 
This  is  achieved  by  specifying  the  forcing  term  Fin  Eq.  (1)  as 
„  d2UB  dUB 

dy  dt  (19) 

where 

dUR_c  dUH 

dt  8  dx 

with  Cg  being  the  disturbance  group  velocity  defined  above. 


3.  NONLINEAR  DEVEIXIP 


vi  I  MicM 


r  OF  GORTLER  VORTICES 


Before  embarking  upon  a  study  of  Gortler/Tollmien-Schlichting  wave 
interaction,  it  is  important  to  study  the  nonlinear  development  of  Gortler 
vortices  to  demonstrate  that  our  numerical  model  adequately  describes  the 
phenomenon  of  interest.  This  we  attempt  to  achieve  by  performing  calcula¬ 
tions  for  the  conditions  of  the  experiment  of  Swearingen  and  Blackwelder 
[14].  The  experiment  of  Swearingen  and  Blackwelder  was  performed  in  a 
low  speed  wind  tunnel  specially  designed  to  have  a  concave  wall  in  order  to 
study  Gortler  instability.  The  radius  of  curvature  of  the  wall  was  r  =  3.2  m. 
Our  calculations  start  at  x  =  45  cm  which,  for  a  free  stream  speed  of  5  m/s, 
gives  an  R  -  -jRex  =  395  and  a  Gortler  number  G  =  7.5.  The  initial 
amplitude  for  the  Gortler  vortex  is  taken  to  be  5%.  The  effective  Gortler 
number  increases  as  the  flow  is  allowed  to  develop  according  to  Eq.  (19). 
The  convection  velocity  Cg  is  taken  to  be  the  disturbance  group  velocity 
which  is  also  used  to  transform  our  temporal  results  to  spatial  results  for 
comparison  with  the  experimental  findings. 
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The  v-w  velocity  vector  plots  are  presented  in  Figure  5  at  computa¬ 
tional  times  of  181  and  880,  the  latter  time  being  quite  far  into  the  nonlinear 
regime.  The  two  plots  are  similar  and  indicate  the  counter-rotating  struc¬ 
ture  of  the  Gortler  vortex  except  that  the  disturbance  intensity  increases 
and  the  vortices  reach  out  further  away  from  the  wall  at  larger  time.  Two 
regions  in  these  figures  may  be  noted:  one  where  the  velocity  vectors  are 
pointed  upwards  away  from  the  wall  and  the  other  where  they  are  pointed 
downwards  to  the  wall.  The  former  will  result  in  a  low  streamwise  velocity 
region  while  the  latter  will  result  in  a  high  velocity  region  as  discussed 
below.  These  regions  may  be  identified  as  "peak"  and  "valley"  following  the 
nomenclature  of  Klebanoff  et  al.  [1]  for  flat  plate  boundary  layer  transition. 

The  distribution  of  the  disturbance  u  velocity  at  various  times 
(distances)  is  plotted  in  Figure  6.  The  solid  lines  represent  positive  velocity 
and  the  dotted  lines  negative  velocity.  Initially  the  two  regions  of  positive 
and  negative  disturbance  velocity  look  similar.  As  the  Gortler  vortex 
develops,  the  region  of  positive  disturbance  velocity  gets  flattened  and 
squashed  near  the  wall  while  the  region  of  negative  velocity  gets  elongated 
in  a  direction  normal  to  the  wall  and  as  a  result  will  influence  the  flow 
further  away  from  the  wall  as  we  note  in  the  next  figure. 

Figure  7  contains  a  comparison  of  the  constant  u  (total)  velocity 
contours  at  various  streamwise  distances.  The  computed  results  are 
presented  in  Figure  7(a)  which  show  the  development  of  mushroom  like 
structure.  The  "mushrooms"  are  caused  by  the  ejection  of  the  low  velocity 
fluid,  away  from  the  surface,  in  the  "peak"  region  noted  above.  Two  peaks 
then  contain  a  "valley"  in  between  as  evident  from  the  figure.  This  valley  is 
caused  by  the  movement  of  high  velocity  fluid  towards  the  wall.  The  exper¬ 
imental  results  for  the  constant  u  velocity  contours  are  given  in  Figure  7(b) 
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and  also  show  the  development  of  mushroom  like  structure.  These  mush¬ 
rooms,  which  are  essentially  the  nonlinear  manifestation  of  the  Gortler 
vortex,  extend  out  of  the  mean  Blasius  boundary  layer  both  in  the  experi¬ 
ment  and  the  present  computations. 

The  wall-shear  stress  along  the  peaks  and  valleys  develop  differently 
as  shown  in  Figure  8  where  comparison  is  made  again  with  the  experi¬ 
ment  of  Swearingen  and  Blackwelder.  In  the  nonlinear  regime,  the  wall 
shear  stress  continues  to  decrease  in  the  "peak”  region  due  to  inflected 
velocity  profiles  while  it  increases  along  the  valley  due  to  the  presence  of 
high  velocity  fluid  near  the  wall.  The  computed  results  would  be  valid  only 
up  to  about  x  =  90  cm,  since  three  dimensional  effects  and  unsteadiness 
caused  due  to  secondary  instability  must  begin  to  appear  which  are  not 
accounted  for  in  the  present  computation.  These  computations  do  show 
that  the  early  nonlinear  stage  is  properly  captured.  A  full  3-D  calculation 
will  be  able  to  capture  the  secondary  instability  stage  and  beyond.  The  3-D 
version  of  the  code  has  been  used  in  the  Gortler/TS  wave-interaction  simu¬ 
lation  and  will  be  discussed  below. 

Due  to  the  inflected  velocity  profiles  in  the  "peak"  region,  the  bound¬ 
ary  layer  considerably  thickens  while  it  thins  in  the  "valley"  region  due  to 
fuller  profiles.  This  is  depicted  in  Figure  9  where  the  computed  boundary 
layer  displacement  thickness  is  compared  with  the  experiment.  Again,  a 
good  agreement  with  the  experiment  is  found  up  to  about  x  =  90  cm. 

The  distribution  of  the  computed  rms  amplitude  of  velocity  distribu¬ 
tion  due  to  Gortler  vortex  is  compared  with  the  experiment  in  Figure  10(a). 
Both  the  experiment  and  the  computations  show  that  the  Gortler  vortex 
reaches  a  saturated  state.  At  this  stage  strong  secondary  instabilities 
would  develop  which  can  be  studied  by  using  our  numerical  model.  We 
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expect  to  observe  at  least  two  types  of  breakdown:  one  involving  a  wavy 
vortex  structure  prior  to  transition  while  the  other  involving  a  horse-shoe 
type  structure.  The  conditions  for  the  appearance  of  these  two  structures 
may  be  established  by  performing  numerical  computations.  Figure  10(b) 
gives  the  disturbance  energy  content  in  various  modes.  The  calculation 
shows  that,  even  though  the  higher  harmonics  grow  rapidly,  most  of  the 
energy  remains  in  the  fundamental  mode.  The  emergence  of  the 
superharmonic  ((0,2)  mode)  is  in  agreement  with  the  experimental 
observation  of  Aihara  &  Koyama  [15]. 

4.  GORTLER/TOT I  MFRN-SCHT  JCHTTNG  WAVE  INTERACTIONS 

The  problem  of  excitation  of  a  Tollmien-Schlichting  wave  in  the 
presence  of  a  Gortler  vortex  is  studied  now.  We  take  the  initial  amplitude  of 
the  Gortler  vortex  to  be  10  percent.  The  non-dimensional  spanwise  wave 
number  is  0.4  and  the  Gortler  number  is  taken  as  12.  Here,  we  use  our 
numerical  model  to  investigate  possible  excitation  of  a  two-dimensional 
(a,0)  TS  wave  in  the  flow  field  modulated  by  the  presence  of  Gortler  vortices. 
The  TS  wave  is  an  upper  branch  mode  and  has  R  =  600,  a  =  .2  and  the  linear 
theory  eigenvalue  is  (.0738,  .0002).  The  results  are  presented  in  Figure  11 
where  (0,1)  mode  indicates  the  Gortler  vortex  and  (1,0)  mode  the  TS  wave. 
The  presence  of  Gortler  vortex  causes  the  growth  rate  of  the  two-dimen¬ 
sional  TS  wave  to  increase  above  its  linear  theory  value.  Towards  the  end  of 
the  computational  domain,  the  growth  rate  of  the  TS  wave  is  higher  than 
the  linear  growth  rate  of  the  primary  (Gortler)  instability. 


The  eigenfunctions  of  the  Gortler  (0,1)  mode  are  presented  in  Figure 
12(a)  for  various  simulation  times.  One  characteristic  that  is  immediately 
evident  is  that  the  shape  of  the  eigenfunction  remains  essentially  the  same 


for  quite  high  amplitudes.  This  was  confirmed  by  restarting  the  calcula¬ 
tions  for  low  initial  amplitudes.  This  implies  that  the  shape  function 
assumption  would  be  applicable  in  a  secondary  instability  analysis  for  the 
Gortler  vortex  problem.  At  later  times,  however,  the  shape  beings  to 
change  and  a  plateau  appears  in  the  eigenfunction  distribution  near  the 
edge  of  the  boundary  layer. 

The  eigenfunction  for  the  TS  (1,0)  mode  are  given  in  Figure  12(b). 
Initially,  only  the  inner  peak  increases  in  amplitude  while  it  moves  slightly 
outwards  from  the  wall.  However,  at  later  times  the  inner  peak  moves 
towards  the  wall  while  increasing  in  amplitude.  After  some  time,  the  outer 
peak  also  increases  in  amplitude  and  becomes  fairly  broad.  The  appear¬ 
ance  of  this  our  structure  is  perhaps  related  to  the  broadening  of  Gortler 
vortex  structure  noted  in  Figure  12(a).  Due  to  the  presence  of  the  Gortler 
vortex  the  mean  flow  is  modulated  and  develops  "wake-like"  profiles  which 
would  then  be  subject  to  instabilities  found  in  laminar  wakes. 

Some  additional  calculations  were  performed  for  a  Gortler  number  of 
14,  p  =  0.5  and  the  initial  Gortler  vortex  amplitude  of  10  percent.  The  TS 
wave  is  again  a  =  .2;  however,  this  wave  decays  linearly  for  R  =  950  which  is 
used  for  the  results  presented  in  Figure  13.  Again  the  (1,0)  mode  grows  in 
the  presence  of  the  Gortler  vortex.  Also  shown  in  the  figure  are  the 
computed  growth  history  of  the  oblique  (1,1)  mode.  Some  numerical  tests 
indicated  that  this  mode  does  not  get  excited  in  the  absence  of  (1,0)  mode. 
However,  this  mode  now  grows  due  to  nonlinear  interaction  of  (0,1)  and 
(1,0)  modes. 

Some  calculations  were  also  performed  at  the  TS  wave  number  a 
which  corresponded  to  the  lower  branch  mode.  Though  more  tests  need  to 
be  performed,  it  appears  that  the  upper  branch  TS  modes  are  excited  while 
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the  lower  branch  modes  are  relatively  insensitive  to  the  presence  of  the 
Gortler  vortex.  Although  the  TS  waves  are  excited  in  the  presence  of  the 
Gortler  vortex  the  amplitude  of  the  Gortler  vortex  required  for  the  interac¬ 
tion  is  quite  high  and  is  at  least  of  the  order  of  10  percent.  Therefore,  the 
Gortler  vortices  will  dominate  the  transition  process  and  the  role  of  the  TS 
waves  will  be  only  in  the  secondary  instability  stage.  However,  here  we 
have  only  considered  the  case  where  Gortler  vortices  are  finite-amplitude 
but  TS  waves  are  small.  The  amplitude  required  for  TS  waves  to  breakdown 
is  much  lower  than  noted  above  for  Gortler  vortices.  Thus  the  situation  will 
be  quite  different  when  TS  waves  are  of  finite  amplitude  or  when  both  are  of 
finite  amplitude.  These  regimes  of  interaction  may  turn  out  to  be  much 
more  dangerous.  It  may  also  be  noted  that  the  presence  of  Gortler  vortices 
may  provide  the  spanwise  scales  for  internalization  of  free  stream 
disturbances  thus  enhancing  the  receptivity  process. 

5.  NONLINEAR  DEVELOPMENT  OF  A  CROSSFLOW  VORTEX 

The  above  methodology  used  for  nonlinear  simulation  of  Gortler 
vortex  problem  equally  applies  to  flows  with  three-dimensional  basic  states, 
such  as  swept-wing  flows  and  flow  due  to  a  rotating  disk.  The  appropriate 
basic  flow  and  initial  conditions,  however,  need  to  be  specified.  The  objec¬ 
tive  here  is  to  carry  out  such  a  calculation  for  a  swept-wing  configuration 
and  see  if  there  are  any  differences  between  the  nonlinear  development  of  a 
crossflow  vortex  in  rotating  disk  and  swept  wing  flow.  The  former  is 
generally  used  as  a  prototype  for  studying  the  instabilities  in  the  swept  wing 
flows.  The  secondary  instability  analysis  of  Reed  [27]  for  a  swept  wing 
suggests  that  there  is  a  resonance  set  up  between  the  crossflow  vortex  and 
its  superharmonic  while  the  computational  results  of  Malik  [28]  for  a  rotat- 
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ing  disk  do  not  suggest  any  such  resonance.  It  is  therefore  important  to 
investigate  any  differences  that  might  exist  between  the  two  flows. 


5.1  Rotating  Disk  Boundary  Laver 

We  first  review  the  results  which  have  been  obtained  for  rotating  disk 
flow.  Just  like  the  swept-wing  flow,  this  boundary  layer  is  amendable  to 
inflectional  instability  which  manifests  itself  in  the  form  of  stationary  co¬ 
rotating  vortices  known  as  crossflow  vortices  since  they  are  caused  by  the 
instability  of  the  crossflow  velocity  profiles.  An  exact  steady  solution  to  the 
Navier-Stokes  equations  exists  for  rotating  disk  flow  and  this  is  one  of  the 
reasons  why  this  flow  has  quite  often  been  used  as  a  prototype  for  studying 
instability  mechanism  in  the  swept  wing  flow.  The  mean  flow  for  the 
rotating  disk  is  governed  by  the  ordinary  differential  (von  Karman  [35]) 
equations: 


F2  -(G  +  l)2  +  F'H-F"  =  0, 


2F(G  +  1)  +  G'H  -  G"  =  0, 
n'  +  HH'-H"  =  0, 

2F  +  H'  -  0, 


(20) 

(21) 

(22) 

(23) 


where  F,  G,  H  and  n  are  related  to  the  velocity  components  and  pressure 
according  to 


Vr  =  rGFiy),  Ve  =  rGG(y), 
Vy=4vGH{y\  P  =  pQvFKy). 


(24) 


where  Q  is  the  angular  speed  of  the  disk.  The  above  equations  are  solved 


subject  to  the  boundary  conditions 
F  =  0,  G  =  0,  H  =  0,  y  =  0 
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(25) 


F  =  0,  G  =  -l  y-»°°  (26) 

A  nonlinear  calculation  was  initiated  by  imposing  a  disturbance 
upon  the  three-dimensional  basic  state  computed  above.  The  disturbance 
shape  was  taken  to  be  that  given  by  the  linear  eigenfunction  which  was 
assigned  an  amplitude  of  10  percent  with  respect  to  the  mean  flow.  The 
wave  numbers  in  the  radial  and  azimuthal  directions  are  a  =  .3751  and 
/3  =  .0777  which  may  be  labelled  as  (1,1)  mode.  The  results  for  the  energies 
of  (1,1),  (2,2),  (3,3),  (4,4)  and  (5,5)  modes  are  presented  in  Figure  14.  The 
solid  curve  represents  the  linear  theory  result  for  the  (1,1)  mode.  Spectral 
broadening  takes  place  and  the  super-harmonics  gain  energy  as  the  growth 
of  the  fundamental  slows  down.  The  energy  cascading  takes  place  through 
(2,2),  (3,3),  (4,4),  etc.  modes.  All  the  modes  tend  to  saturate  and  the 
harmonics  do  not  overtake  the  fundamental  mode. 

The  azimuthal  velocity  profiles  at  fixed  radius  but  at  various 
azimuthal  locations  within  the  computational  domain  are  presented  in 
Figure  15.  The  nondimensional  time  of  the  simulation  is  170.  While  the 
von  Karman  profile  is  full  with  no  inflection  point,  the  computed  flow 
develops  "wake-like"  profiles  at  y  around  2.5.  These  profiles  are  amenable 
to  secondary  instabilities.  However,  no  attempt  was  made  here  to  capture 
these  instabilities. 

Hot-wire  traces  at  various  radii  from  the  experiment  of  Wilkinson 
and  Malik  [22]  are  presented  in  Figure  16(a).  The  hot-wire  is  located  at 
y  =  1.8  and  is  oriented  such  that  it  essentially  senses  the  perturbations  in 
the  azimuthal  velocity  component.  The  traces  develop  a  "kink"  which  may 
be  attributed  to  the  superharmonic.  The  computational  traces  at  various 
radii  are  presented  in  Figure  16(b).  The  nondimensional  time  of  the  simu- 


19 


lation  is  again  170.  The  computational  trace  appears  to  develop  the  same 
type  of  "kink"  as  in  the  experiment.  The  amplitude  of  this  "kink",  appar¬ 
ently  caused  by  the  superharmonic,  is  small  compared  to  the  fundamental 
both  in  the  computation  and  the  experiment.  The  superharmonic  appears 
due  to  the  self-interaction  of  the  crossflow  vortex. 


.2  Falkner-Skan-Cooke  Boundary  Laver 

Basic  insight  into  instability  of  flow  past  a  swept  wing  may  be 
achieved  by  considering  Falkner-Skan-Cooke  (FSC)  flow.  To  this  end,  let  us 
consider  boundary  layer  flow  over  an  infinite  swept  wedge  at  zero 
incidence.  The  governing  equations  may  be  written  as  (see  Rosenhead  [36]) 


f'"  +  ff"  +  P  (l-f")  =  0 

(27) 

g"+fg'  =  Q 

(28) 

where  ul  ue  -  f\vp,  wl  we  =  g{ri)  (ue  and  we  being  the  boundary  layer  edge 
velocities  in  the  chordwise  x  and  spanwise  z  directions  (see  Figure  17) 
respectively),  and  p  is  the  Hartree  pressure  gradient  parameter.  The 
similarity  parameter  r\  is  defined  as 


V  = 


u. 


(2- fi)vx 


(29) 


The  above  equations  are  to  be  solved  subject  to  the  boundary 
conditions: 


f  =  f'  =  g  =  0  ;  7]  =  0  (30) 

f'-*1  S  -*  1  ;  n-»°°  (31) 

With  reference  to  Figure  17,  the  streamwise  and  crossflow  velocity 
profiles  can  be  constructed  for  FSC  flow  as 

U(y)  =  f'(y)  cos2  6  +  g(y )  sin2  6  (32) 
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WXy)  =  [-/"'O') +  £(:y)]sin0cos0  (33) 

where  both  the  streamwise  and  crossflow  velocities,  U  and  W  in  xs  and  zs 


directions  respectively,  have  been  scaled  with  Ve  =  ■\fuf+wf  and  the  external 
streamline  angle  6  is  given  as 
6  =  tan~'( 


\U'J 


(34) 


The  functions  U  and  W  for  (i  =  0.6  are  ploltted  in  Figure  18  where  y  is 
the  nondimensional  distance  normal  to  the  surface  and  has  been  scaled 
with  l  =  ^  vx  /  ue . 

In  Figure  17,  xu  is  along  the  axis  of  the  crossflow  vortex  which  is 
oriented  at  some  small  angle  <p  with  the  free  stream  direction.  The  vortex 
coordinate  system  (x0,y,zu)  is  used  for  the  nonlinear  simulation  of  the 
crossflow  vortex.  The  analysis  begins  by  imposing  a  disturbance  of  the 
form 

u(xu,y,zv,t)  =  u(y)eilPA~oM} 


on  the  FSC  flow.  The  wave  front  is  then  aligned  with  xu  direction.  Here  we 
consider  a  stationary  disturbance  so  that  Real  (a>)  =  0  and  Im(co) 
corresponds  to  the  temporal  growth  rate  of  the  vortex.  The  calculations  are 
performed  for  the  following  conditions: 


R  =  Vtl/  v  =  800 

and 

A  =  0.48,  <p  =  -4.135°. 


The  above  combination  of  parameters  results  in  a  stationary  vortex 
(cor  =  0)  disturbance  according  to  linear  theory.  The  linear  theory  growth 
rate  is  .0065.  The  crossflow  Reynolds  number  is  244. 

Linear  eigenfunctions  were  first  obtained  to  provide  initial  conditions 
for  the  nonlinear  simulations.  The  eigenfunctions  uyvyw  are  plotted  in 
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Figure  19.  These  can  be  compared  with  the  results  presented  in  Figure  1 
for  Gortler  vortex  which  is  also  a  stationary  streamwise  vortex.  We  note 
that  the  eigenfunctions  are  complex  for  the  crossflow  vortex;  however,  just 
like  the  Gortler  vortex  the  magnitudes  of  v,w  components  are  much 
smaller  than  the  u  component  along  the  axis  of  the  stationary  crossflow 
vortex. 

These  eigenfunctions  were  included  in  the  initial  conditions  for  the 
numerical  study  of  the  nonlinear  evolution  of  the  crossflow  vortex.  The  ini¬ 
tial  amplitude  was  taken  to  be  5  percent.  The  ( v,w )  vector  plots  are  pre¬ 
sented  in  Figure  20.  As  opposed  to  the  vector  plots  given  in  Figure  5  for  the 
Gortler  vortex,  it  can  be  seen  that  crossflow  vortices  are  corotating  vor+ices 
with  one  vortex  per  wave  length.  The  disturbance  intensity  increases  as  the 
vortices  develop  nonlinearly  (compare  Figures  20(a)  and  (c)). 

The  evolution  of  the  energies  of  various  modes  is  given  in  Figure  21. 
Just  like  in  the  case  of  a  Gortler  vortex,  higher  harmonics  grow  and  they 
all  tend  to  saturate.  The  energy  in  the  fundamental  (0,1)  mode  dominates. 
This  is  quite  similar  to  what  was  observed  above  for  a  crossflow  vortex  in 
rotating  disk  flow. 

Saric  &  Yeates  [26]  observed  that,  in  addition  to  the  fundamental 
wavelength,  an  additional  peak  in  the  hot-wire  trace  appears  at  some 
height  from  the  wall  with  half  the  wavelength.  This  latter  peak  had  an 
amplitude  about  3  times  that  of  the  fundamental.  The  surface  flow  visual¬ 
izations  only  showed  streaks  at  the  wavelength  given  by  the  fundamental. 
Reed  [27]  attempted  to  explain  the  anomaly  found  in  the  experiment  by 
using  secondary  instability  theory  based  upon  Floquet  analysis.  Thus  the 
periodic  basic  flow  was  constructed  by  superposing  on  the  three-dimen¬ 
sional  boundary  layer  the  eigenfunctions  for  the  crossflow  vortex  of 
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wavenumber  pu.  The  amplitude  assigned  to  this  disturbance  was  7  percent 
in  order  to  simulate  the  experimental  conditions.  Linear  instability  analy¬ 
sis  of  this  basic  flow  revealed  that  a  disturbance  of  wavenumber  2/?  is 
amplified  while  this  disturbance  is  stable  from  a  linear  analysis  of  the 
boundary  layer  flow  alone.  However,  the  growth  rate  for  the  superhar¬ 
monic  2f}u  was  found  to  be  almost  the  same  as  for  the  fundamental  (Sv  mode 
and,  thus,  one  of  the  basic  assumptions  of  the  secondary  instability  theory 
was  violated.  In  any  case,  the  eigenfunctions  obtained  from  the  secondary 
instability  analysis  were  assigned  an  amplitude  of  21  percent  (3  times  the 
fundamental)  and  the  streamlines  for  the  resulting  flow  were  constructed 
which  showed  two  vortices  per  wavelength  of  the  fundamental  instability. 
The  results  appeared  to  explain  the  observations  of  Saric  and  Yeates  [26], 
However,  we  note  that  in  our  nonlinear  calculations  presented  in  Figure  21, 
the  superharmonic  does  not  achieve  amplitude  higher  than  that  of  the 
fundamental,  let  alone  three  times  the  amplitude  of  the  fundamental. 

It  may  well  be  that  the  superharmonic  has  an  eigenfunction  struc¬ 
ture  which  has  local  peak  somewhere  in  the  boundary  layer  at  an  ampli¬ 
tude  of  about  3  times  that  of  the  fundamental  as  observed  in  the  experiment 
of  Saric  and  Yeates.  Therefore,  we  plot  the  eigenfunctions  for  the  funda¬ 
mental  and  the  superharmonic  at  four  different  times  in  Figure  22.  The  u 
amplitude  of  the  fundamental  and  the  superharmonic  increases  with  time 
but  no  where  in  the  boundary  layer  the  latter  attains  amplitude  higher  than 
that  of  the  former. 

The  fundamental  u  perturbation  eigenfunction  develops  a  "kink" 
which  at  later  times  develops  into  a  double  peak.  Somewhat  similar  struc¬ 
ture  appears  in  the  swept  wing  experiments  of  Dagenhart  et  al.  [25]  but  is 
unsteady.  At  this  time,  we  have  not  attempted  to  find  the  frequency  content 
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of  the  computed  signal.  Such  a  comparison  will  be  meaningful  when  the 
calculations  are  performed  for  the  exact  conditions  of  the  experiment  of 
Dagenhart  et  al.  We  also  note  that  initially  the  shape  of  the  fundamental 
eigenfunction  is  preserved  even  though  the  amplitude  is  relatively  lavge. 
This  lends  some  support  to  the  shape  function  assumption  generally  used 
in  the  secondary  instability  analysis. 

The  computational  traces  of  u  at  various  heights  in  the  boundary 
layer  and  at  three  different  times  are  presented  in  Fig.  23.  The  u{.zu)  velocity 
develops  a  "kink"  similar  to  that  found  for  the  rotating  disk  case  due  to  the 
development  of  the  superharmonic.  This  kink  is  quite  pronounced  at  a 
certain  height  from  the  surface.  A  hot-wire  traversed  in  the  spanwise 
direction  at  this  particular  height  will  indicate  the  presence  of  a  distur¬ 
bance  with  half  the  wavelength  of  the  fundamental;  however,  as  noted 
above  the  amplitude  of  the  superharmonic  is  less  than  the  fundamental. 
The  observation  of  Saric  &  Yeates  that  the  superharmonic  attains  three 
time  the  amplitude  of  the  fundamental  may  be  the  result  of  certain  condi¬ 
tions  unique  to  the  experiment,  particularly  when  the  experiments  of 
Muller  &  Bippes  [24]  and  Dagenhart  et  al.  [25]  fail  to  observe  the 
phenomenon  noted  by  Saric  &  Yeates.  It  may  be  that  the  nonparallel/ 
nonlinear  calculations  for  the  exact  conditions  of  the  Saric-Yeates  experi¬ 
ment  will  explain  the  observed  results.  Nonparallel  effects  are  accounted 
for  neither  in  the  secondary  instability  analysis  of  Reed  [27]  nor  in  the 
present  temporal  numerical  simulations.  In  any  case,  we  conclude  that 
the  anomaly  found  in  the  Saric-Yeates  experiment,  namely  that  the  super¬ 
harmonic  attains  an  amplitude  three  times  that  of  fundamental,  remains 
unexplained. 
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The  velocity  profiles  parallel  and  at  90°  to  the  crossflow  vor tex  axis  at 
various  spanwise  locations  are  presented  in  Figure  24.  Again,  inflected 
"wake-like"  profiles  develop  in  the  direction  of  the  vortex  axis  (which  is  at 
an  angle  of  -4.135°  with  respect  to  the  inviscid  streamline)  due  to  nonlinear 
distortion  of  the  crossflow  vortex  as  also  found  for  the  rotating  disk.  These 
profiles  will  be  subject  to  inviscid  secondary  instability.  From  the  results 
presented  in  this  section,  it  becomes  quite  clear  that  the  nonlinear  develop¬ 
ment  of  a  crossflow  vortex  in  FSC  boundary  layer  is  quite  similar  to  what 
was  found  for  a  rotating  disk. 

5.3  Secondary  Instability  in  FSC  Boundary  Laver 

The  computed  Falkner-Skan-Cooke  flow  field  in  the  presence  of  a 
crossflow  vortex  has  been  presented  above.  It  is  of  interest  to  investigate  the 
type  of  secondary  instabilities  this  flow  may  be  subjected  to.  The  present 
numerical  simulations  have  been  used  to  investigate  one  such  instability. 
The  calculations  have  been  performed  for  the  parameter  values  used  above 
in  the  nonlinear  simulation  of  the  crossflow  vortex.  However,  included  in 
the  initial  conditions  is  a  mode  with  a v  =  0.48  and  /?,,  =  0.48.  This  unsteady 
mode  is  highly  damped  according  to  the  linear  theory'.  The  growth  rate  of 
this  mode  computed  from  our  numerical  simulations  is  plotted  in  Figure 
25.  Initially,  the  mode  is  damped,  in  agreement  with  the  linear  stability 
theory.  However,  after  about  t  =  100,  this  mode  begins  to  grow.  The  growth 
rate  distribution  is  oscillatory  because  of  the  interactions  with  other  modes, 
but  eventually  the  growth  rate  of  this  mode  is  much  larger  than  the  linear 
growth  rate  of  the  crossflow  vortex  ((0,1)  mode)  indicated  by  circles  in  the 
figure.  The  computed  growth  rate  of  the  (0,1)  mode  indicated  by  a  dashed 
line  is  also  shown  and  it  approaches  zero  as  the  crossflow  mode  saturates. 
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The  (a„,/J„)  mode  (denoted  as  (1,1)  mode)  thus  develops  quite  quickly  and  will 
perhaps  provide  at  least  one  mechanism  for  breakdown  of  the  crossflow 
vortex.  The  values  of  au  and  (5V  chosen  here  imply  a  wave  oriented  at  an 
angle  of  about  45°  with  respect  to  the  crossflow  vortex. 

For  rotating  disk  flow,  Kohama  [37]  observed  that  prior  to  transition  a 
new  instability,  in  the  form  of  small  vortices,  appears  on  the  surface  of 
crossflow  vortices.  These  vortices  were  oriented  at  an  angle  of  about  45° 
with  respect  to  the  crossflow  vortex.  The  structure  observed  by  Kohama  is, 
perhaps,  the  result  of  a  localized  instability  which  takes  place  at  some 
distance  away  from  the  wall  and  produces  an  illusion  of  the  "vortices" 
riding  on  the  crossflow  vortex.  Obviously,  the  flow  visualization  technique 
used  can  not  provide  the  true  nature  of  the  instability.  In  any  case,  it 
appears  that  the  secondary  instability  computed  here  may  provide  a  clue  for 
the  explanation  of  the  experiment  of  Kohama.  It,  therefore,  would  be  desir¬ 
able  to  further  investigate  this  phenomenon  using  Floquet  analysis  and 
more  numerical  simulations.  It  may  be  noted  here  that  since  the  growth 
rate  of  the  (av,(5v)  mode  is  about  4-5  times  the  linear  growth  rate  of  the  cross- 
flow  vortex  ((0,/IJ  mode),  Floquet  analysis  can  provide  useful  results. 
Further  investigations  into  the  nature  of  this  secondary  instability  will  be 
made  at  a  later  date. 

6.  CONCLUSIONS 

We  have  used  direct  numerical  simulations  to  investigate  nonlinear 
development  of  Gortler  and  crossflow  vortices  and  to  study  Gortler/TS  wave 
interactions.  Our  results  show  that  when  Gortler  vortex  develops  nonlin- 
early,  higher  harmonics  grow  but  most  of  the  energy  remains  in  the 
fundamental  mode.  The  flow  field  in  the  presence  of  the  nonlinear  Gdrtler 
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vortex  has  been  computed  for  the  conditions  of  the  experiment  of 
Swearingen  &  Blackwelder  [14].  In  agreement  with  the  experiment,  a 
"mushroom"  like  structure  appears  as  the  low  velocity  fluid  is  ejected  up¬ 
wards  between  the  counter-rotating  Gortler  vortices.  In  this  "peak"  region, 
wall  shear  stress  decreases  while  the  displacement  thickness  increases 
relative  to  the  mean  Blasius  boundary  layer.  The  converse  is  true  for  the 
"valley"  region  where  the  high  velocity  fluid  is  pushed  towards  the  wall. 

Numerical  simulations  are  also  used  to  study  the  excitation  of  a  two- 
dimensional  TS  wave  in  the  presence  of  a  finite-amplitude  Gortler  vortex. 
The  results  indicate  that  the  TS  wave  grows  faster  than  its  linear  growth 
rate  and  that  oblique  TS  waves  also  grow  due  to  nonlinear  interactions. 
These  oblique  waves  are  linearly  stable.  The  amplitude  of  the  Gortler  vortex 
at  which  the  interaction  takes  place  are  much  larger  than  that  needed  in 
flows  where  finite-amplitude  TS  waves  constitute  the  primary  instability. 

Finally,  we  have  used  the  computational  approach  to  simulate  the 
nonlinear  development  of  a  crossflow  vortex  in  a  Falkner-Skan-Cooke 
boundary  layer.  It  is  shown  that  nonlinear  development  of  a  crossflow 
vortex  in  FSC  and  rotating  disk  boundary  layer  are  similar.  In  both  flows, 
the  fundamental  saturates  at  sufficiently  large  amplitudes  along  with  the 
harmonics.  The  first  superharmonic  clearly  emerges  in  the  velocity  signal 
but  its  amplitude  remains  smaller  than  the  fundamental.  A  new  type  of 
secondary  instability  is  noted  in  the  numerical  simulations.  This  instabil¬ 
ity  is  found  to  have  growth  rates  which  are  4  to  5  times  higher  than  the 
growth  rate  of  the  linear  crossflow  vortex. 
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Figure  1.  Eigenfunctions  for  a  Gortler  vortex  in  a  Blasius  boundary  layer, 
G=  14,/?  =  .3. 
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Figure  2.  Eigenfunctions  for  an  oblique  TS  wave  in  a  Blasius  boundary 
layer,  R  =  950,  a  =  .103,  P  =  .15. 


32 


Figure  3. 


TIME 


Computed  evolution  of  energies  of  a  two-dimensional  TS  wave 
((a, 0)  mode)  and  the  oblique  subharmonic  ((ct/2,/?  mode).  In  this 
calculation,  R  =  606,  a  =  .2033,  /3  =  .16968.  A  strong  secondary 
subharmonic  instability  develops.  Included  in  the  initial 
conditions  are  two  oblique  primary  subharmonics.  The  solid 
lines  represent  linear  theory  results. 
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Figure  4.  Same  as  for 
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TtMC  —  15.203  OlST  —  48,872 


(a)  Time  =  15.2 
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(c)  Time  =  766 

Figure  6.  Disturbance  u-velocity  in  y-z  plane  for  Gortler  vortex  flow  at 
various  times. 
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Figure  7(a).  Computed  iso-contours  of  the  streamwise  velocity  in  the  cross 
stream  y-z  plane  for  Cortler  vorte  .  flow.  The  contours  are 
plotted  with  velocities  0.1,  0.2,. ..,0.9. 
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Figure  10(b).  Evolution  of  disturbance  energies  of  various  Gortler  modes  in 
a  boundary  layer.  The  fundamental  (0,1)  mode  continues  to 
have  the  dominant  energy. 
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Figure  11.  Evolution  of  energies  of  various  modes  using  direct  numerical 
simulation  and  comparison  with  linear  theory  results.  Due  to 
the  presence  of  Gortler  vortex  ((O,l)mode),  the  TS  wave  ((1,0) 
mode)  grows  faster  than  its  linear  theory  growth  rate. 
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Figure  14.  Computed  evolution  of  energies  of  various  modes  in  the  presence 
of  a  finite-amplitude  crossflow  vortex  ((1,1)  mode)  in  rotating 
disk  flow.  In  this  calculation  R  =  500,  a  =  .3751,  B  =  .0777.  Initial 
amplitude  of  the  (1,1)  mode  is  10%  of  the  mean  flow. 
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Figure  15.  Computed  azimuthal  velocity  pi-ofiles  in  rotating  disk  flow  at  a 
fixed  radial  location  and  within  1/2  fundamental  wavelength. 
Time  t  =  170. 
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Figure  16(a)  Hot-wire  traces  at  various  radii  near  transition  in  rotating 
disk  boundary  layer.  Hot-wire  is  located  at  y  =  1.8  and  it  is 
oriented  such  that  it  senses  the  perturbation  in  azimuthal 
velocity  component. 
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Figure  16(b)  Computed  traces  of  azimuthal  velocity  perturbations  at 
y  =  1.8.  Only  one  fundamental  wa/elength  is  shown.  The 
calculations  are  for  the  conditions  of  Fig.  14.  Time  l  =  170. 
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Figure  17.  Coordinate  system  for  Falkner-Skan-Cooke  boundary  layer. 
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(a)  velocity  component  along  the  axis  of  the  vortex 


Figure  19.  Disturbance  eigenfunctions  for  a  crossflow  vortex  in  FSC 
boundary  layer. 
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(b)  velocity  component  normal  to  the  surface 
Figure  19.  Continued. 
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(c)  velocity  component  in  the  direction  of  the  wavenumber  vector 

Figure  19.  Concluded. 
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Figure  20.  Continued 
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Figure  20.  Concluded. 
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(a)  Tii 

Figure  22.  Disturbance  u  profiles  for 
at  various  times  in  the  pre 
vortex. 
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(d)  Time  =  327 
Figure  22.  Concluded. 


64 


y 

4Sf) 


Figure  23. 
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(b)  Time  =  195 
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(c)  Time  =  355 

Computational  trace  of  the  velocity  component  along  the  axis  of 
the  crossflow  vortex  at  three  different  times  and  at  various 
heights  from  the  surface  for  FSC  boundary  layer.  Two  wave 
period  are  shown. 
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Figure  24.  Concluded 


